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New blazars from the cross-match of recent 
multi-frequency catalogs 
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Abstract Blazars are radio-loud active galactic nuclei 
well known for their non thermal emission spanning a 
wide range of frequencies. The i?oma-BZCAT is, to 
date, the most comprehensive list of these sources. We 
performed the cross-match of several catalogs obtained 
from recent surveys at different frequencies to search 
for new blazars. We cross-matched the I'** Swift XRT 
Point Source catalog with the spectroscopic sample of 
the 9*^ Data Release of the Sloan Digital Sky Survey. 
Then, we performed further cross-matches with the cat¬ 
alogs corresponding to the Faint Images of the Radio 
Sky at Twenty cm survey and to the AllWISE Data re- 
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lease, focusing on sources with infrared colors similar to 
those of confirmed y-ray blazars included in the Second 
Fermi-LAT catalog. As a result, we obtained a prelim¬ 
inary list of objects with all the elements needed for a 
proper blazar classification according to the prescrip¬ 
tions of the Roma-BZCAT. We carefully investigated 
additional properties such as their morphology and the 
slope of their spectral energy distribution in the radio 
domain, the features shown in their optical spectrum, 
and the luminosity in the soft X rays to exclude generic 
active galactic nuclei and focus on authentic blazar-like 
sources. At the end of our screening we obtained a list 
of 15 objects with firmly established blazar properties. 

Keywords X-rays: galaxies - galaxies: active - radia¬ 
tion mechanisms: non-thermal. 


1 Introduction 


Blazars are peculiar radio-loud active galactic nuclei 
(AGNs) emitting over the whole electromagnetic spec¬ 
trum, from the radio band to the MeV-GeV range of 
7 rays and extending towards the TeV range in the 
most extreme cases. Large-amplitude flux and spectral 
varia bility, on short time scales down to a f ew hours 
{i.e., Montagni et al. 200d Guota et al. 2008 ). charac¬ 
terize their emission. In the radio band they appear 
as core-dominated, flat spectrum so urces. The discov¬ 


ery of q-ray emission frorn blazars (ISwanenburg et al 


1978 : Hartman et al. 1999t Ackermann et al. 2011 ) es¬ 


tablished them as a class of extreme particle accel¬ 
erators which generate some of the most energetic 
photons of any known extragalactic cosmic source. 
Within the framework of the AGN Unification Scenario 


(|Urrv and Padovanil 119951) . these properties are inter¬ 
preted as non thermal radiation emitted by blobs of 
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charged particles moving away from a central super- 
massive black hole at velocities close to the speed of 
light and in directions very close to the line of sight. 
In these conditions, the non thermal radiation is rel- 
ativistically amplified and becomes dominant with re¬ 
spect to all other components of thermal origin. The 
blazar class includes both BL Lac objects (BL Lacs) and 
Flat Spectrum Radio Quasars (FSRQs), similar in the 
shape of their spectral energy distribution (SED) but 
different in their optical spectrum. BL Lacs display 
no or very we ak emission lines w i th equivalent widt h 


Wa ^ 5 a ( Stickel et al.lll991t Istocke et al.lll991 ): 


conversely, FSRQs reveal evidence of broad emission 
lines. According to the Unification Scenario, BL Lacs 
and FSRQs are thought to be the beamed counterparts 
of low- (FRI) and high-lum inosity (FRII) radio galaxies 


( Fanaroff and Rilevlll974ll . respectively. 


In recent years the number of all-sky surveys in dif¬ 
ferent energy bands has considerably increased. This 
large amount of data is very helpful for blazar research 
but needs to be properly managed. At this purpose a 
few, well-defined and stringent criteria for blazar classi- 
fication have been defi ne d lead i ng to t he Roma-BZCAT 
( Massaro et al. 2009al f9. 20ld 2011a). These criteria 


are: 1) detection in the radio band down to mJy flux 
densities at 1.4 GHz; 2) for FSRQs, a radio spectral 
index ar lower than a* = 0.5 (F(:/) oc v~'^) measured 
between 1.4 GHz and 5 GHz: this condition is not re¬ 
quired for BL Lacs, although most of them have flat 
spectra; 3) compact radio morphology or, when ex¬ 
tended, with a dominant core and a one-sided jet; 4) 
optical identification and analysis of the optical spec¬ 
trum to establish the type of blazar; 5) isotropic X-ray 
luminosity Lx > L\ = 10"^^ erg s“^ from a point-like 
source. The 4*^ version of the Roma-BZGAT, available 
onlin^, included 3,149 objects: 1,221 BL Lacs (BZBs in 
the catalog, -38.8%), 1,707 FSRQs (BZQs, -54.2%), 
and 221 blazars with uncertain classification (BZUs, 
—7.0%). 270 sources were included among the BL Lacs 
despite the lack of an optical spectrum and were there- 


campaigns (e.Q 

.. Landoni et al. 

2013; 

Pae-ffi et al. 

2014; 

Massaro et al. 

20141 have been carried out to confirm 


blazars (ID’Abrusco et al.ll2012fl led to the parametriza- 
tion of the region that they occupy in the color space 
of the Wide - field Infrared Survey Explorer (WISE, 


Wright et al.l 1201011. named the WISE blazar locus 


(jMassaro et al.N2012ll . This parametrization w as used 


to develop a method (jP’Abrusco et al.l 1201311 for the 


association of unidentified y-ray sources detected by 
Fermi-LPLL I Atwood et al.l [200911 and included in the 
Second Fermi-LAT catalog f2FGL. IXolan et"^ 2012 1. 
According to this method, the position of a generic 
IR source with respect to the locus can be considered 
a useful indicator of its blazar nature. Recently, the 
parametrization has been slightly modified to take into 
account the results of the AllWISE Data Releas^ lead¬ 
ing to the compilat ion of an all-sky catal og of candi¬ 
date y-ray blazars (ID’Abrusco et al.N2014l l. This cat¬ 
alog includes the WISE blazar-like radio-loud sources 
(WIBRaLS) detected in all the four WISE filters. 

In this work we provide a list of sources with firmly 
established blazar properties after a critical review of 
a preliminary list of objects obtained from the cross¬ 
match of recent multi-frequency surveys. The paper 
is organised as follows: the guidelines adopted in the 
cross-match of the survey catalogs used in this work 
are illustrated in Section!^ The validation of a sample 
of blazars among a preliminary list of candidates is pre¬ 
sented in Section [H Finally, a discussion of our results 
is reported in Section SI 

The WISE magnitudes in the [3.4], [4.6], [12], and 
[22] ^m nominal filters are expressed in the Vega 
system; the magnitude values of three WISE filters 
(namely those corresponding to [3.4], [4.6], and [12] ^m) 
and of the colors [3.4]-[4.6], [4.6]-[12], and [12]-[22] de¬ 
rived using these magnitudes, are corrected for Galac¬ 
tic extintio n acco rding to the extinction law presented 
by Drainel ( 2003l l. Unless stated otherwise, we used 
c.g.s. units throughout and assumed a A cold dark mat¬ 
ter (AGDM) cosmology wit h Hn = 72 km s~^ M pc~^, 
Dm = 0.26, and Hy = 0.74 ( Dunklev et al. 2009). 


2 Cross-match of the used survey catalogs 


the BL Lac classification of these sources; additional We started our analysis taking into a ccount, amon; 
campaigns are ongoing. 

New, complementary approaches to the blazar clas¬ 
sification come from the evaluation of the ir color in¬ 
dices in the infrared ( Massaro et al. 2011bll as well as 


in the optical-UV bands (jMassaro et al.l l2012h . The 


discovery of peculiar infrared colors of y-ray emitting 


the cr iteria for blazar’s classification (jMassaro et al.l 
2009bll . those related to the point-like X-ray emission 
and to the possibility of analysing the optical spec - 
trum. Therefore, we used TOPGAtH (iTavloij l2(105ll 
to cross-mat ch the 1^* Swift X RT Point Source cata¬ 
log llSXPS. lEvans et al.ll2014l) with the spectroscopic 


^http://www.asdc.asi.it/bzcat4/ 


^http: //wise2. ipac.caltech.edu/docs/release/allwise/expsup 
^http: //www. star.bris.ac.uk/'^mbt/topcat/ 

































































































3 


Table 1 The list of validated blazars 


name 

dox’^ 

Rx'" 

doR° 

doi"^ 

S1.4' 



r 

z 

Fx'^ 


classification 

J003931.58-111102.4 

1.2 

4.2 

0.41 

0.24 

296 

237 

-HO.18 

19.0 

0.553 

6.30 

75.1 

FSRQ 

J082438.99+405707.7 

2.1 

4.7 

0.51 

0.09 

175 

138 

-HO.19 

17.0 

0.612 

0.10 

1.5 

FSRQ 

J082753.69+521758.3 

2.6 

3.8 

0.68 

0.14 

181 

292 

-0.38 

18.9 

0.338 

2.50 

9.2 

FSRQ 

J095906.96+050958.9 

1.5 

4.4 

0.42 

0.09 

47 

73 

-0.35 

18.0 

0.996 

0.38 

19.5 

FSRQ 

J100033.84+132410.8 

1.7 

3.9 

0.20 

0.13 

53 

34 

-HO.36 

16.7 

1.355 

0.74 

82.2 

FSRQ 

J103045.22-^255522.1 

3.1 

5.1 

0.18 

0.04 

49 

50 

-0.02 

17.1 

0.692 

0.24 

5.0 

FSRQ 

J104031.62+061721.7 

1.1 

3.9 

0.12 

0.13 

39 

49 

-0.18 

19.9 

0.743? 

0.66 

15.9 

BL Lac 

J110838.98+255613.2 

0.3 

3.8 

0.20 

0.31 

68 

73 

-0.06 

17.5 

0.732 

0.80 

19.0 

FSRQ 

J133631.44+031423.5 

5.4 

4.3 

0.32 

0.75 

62 

111 

-0.47 

18.9 

1.303 

0.12 

12.1 

FSRQ 

J141238.66-^484447.1 

1.6 

3.8 

0.16 

0.13 

58 

74 

-0.20 

18.6 

0.906 

0.39 

15.8 

FSRQ 

J142114.05-^282452.8 

3.4 

4.1 

0.68 

0.03 

49 

40 

-HO.16 

17.8 

0.538 

0.46 

5.1 

FSRQ 

J153458.41+575625.6 

2.5 

4.6 

0.08 

0.29 

7 

17 

-0.71 

19.0 

1.129 

0.17 

11.9 

FSRQ 

J161541.21-^471111.7 

2.9 

4.6 

0.23 

0.09 

98 

175 

-0.46 

17.2 

0.199 

1.28 

1.4 

BL Lac 

J162805.20+252636.8 

1.8 

4.7 

0.14 

0.08 

100 

59 

+0.42 

18.4 

0.995 

0.31 

15.9 

FSRQ 

J170634.12+361508.0 

4.9 

4.4 

0.37 

0.14 

19 

27 

-0.28 

17.9 

0.917 

1.06 

44.2 

FSRQ 


•^Angular separation between the SDSS and the ISXPS sources in arcseconds 
^90% error radius of the ISXPS source in arcseconds 

^Angular separation between the SDSS source and the corresponding FIRST associated counterpart in arcseconds 
‘^Angular separation between the SDSS source and the corresponding WISE associated counterpart in arcseconds 
^Radio flux density at 1.4 GHz in mjy 
•f Radio flux density at 5 GHz in mJy 
S'Radio spectral index cxr 

^Unabsorbed X-ray flux (0.3-10 keV) in 10“^^ erg cm“^ s“^ 

^X-ray luminosity (0.3-10 keV) in 10“^^ erg s“^ 


sample of the 9*^ Data Release of the Sloan Digital 
Sky Survey (SDSS-DR9, Ahn et al. 2012). The ISXPS 
catalog resulted from the analysis of data collected by 
the Swift X-Ray Telescope (XRT, ISurrows et al. 200511 
across 8 years of operations since its launch, up to 
October 2012. It included 151,524 soft X-ray point¬ 
like sources detected in the 0.3-10 keV band with an 


appro ximately uniform sky distribution. lEyans et al. 


1 2 OI 4 I 1 carried out seyeral checks to filter out spuri¬ 
ous or extended sources and to express the reliability 
of each detected source. As a result, a quality flag 
which can assume three integer yalues (0, 1, and 2) 
to indicate a good, reasonable, or poor detection, re- 
spectiyely, was assigned to each source. In our cross¬ 
match we considered sources from a cleaned yersion of 
the ISXPS including only detections flagged as good 
or reasonable, made of 98,762 objects. We spatially 
cross-matched this list with the SDSS-DR9 suryey cat¬ 
alog searching for optical counterparts, with ayailable 
spectra, within the error circle (at 90% confidence) of 
each X-ray source. At this purpose, we added to this 
90% error radius a positional uncertainty of 2" to take 
into account the typical seeing of the SDSS. The list 
of 5012 SDSS sources that we obtained was then cross- 


matc hed with the WIBRaLS catalog (ID’Abrusco et al 


2 OI 4 I 1 . We use d an oyerall matching r adius of 3.3" as 


determined by D’Abrusco et al. ( 2013h in their search 
for the optimal yalue of the spatial association between 


blazars in the Roma-BZGAT and WISE sources, and 
obtained 214 objects. All the WIBRaLS haye a radio 
counterpart. Howeyer, we kept for further inyestigation 
only objects with a radio counterpart listed in the Faint 
Images of the Radio sky at Twenty cm Suryey (FIRST, 


Becker et al.l Il995t IWhite et al.l 1199711 since the angu¬ 


lar resolution of this suryey allows a better eyaluation 
of the radio morphology. In the cross-match between 
our list of 214 SPSS sources and the l atest yersion of 
the FIRST catalog ( Becker et al. 20I2l l we adopted a 
matching radius of 4" to take into account the posi¬ 
tional uncertainty on the FIRST source (2") in addi¬ 
tion to the SDSS one. We obtained 201 objects includ¬ 
ing 152 blazars already present in the 4‘^ yersion of 
the i?omo-BZCAT: consequently, we achieyed a prelim¬ 
inary list of 49 blazar candidates whose properties were 
inyestigated in greater detail. 


3 Validation of blazar candidates 

The two main aspects that we took into account in 
the examination of the radio emission properties of our 
candidates were the morphology and the slope of their 
SED. Blazars are characterized by a core-dominated ra¬ 
dio morphology; as an alternatiye, the additional pres¬ 
ence of a one-sided jet can be accepted. The eyi- 
dence of extended structures, and in particular those 
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shaped in the form of two lobes emerging from the 
central core following an FRI or an FRII morphology 
( Fanaroff and Rilev 1974ll is the clear footprint of a ra¬ 
dio galaxy instead of a blazar. Other irregular extended 
structures, adequately resolved by the FIRST survey, 
were not accepted either. As a result, we excluded from 
our preliminary list 11 radio galaxies, mainly character¬ 
ized by a FRII morphology. 

Then we estimated the slope in the SED at radio 
frequencies of the remaining 38 objects finding a steep 
spectrum for 13 of them. In our list we also found 10 
objects with no measurements of flux density S at radio 
frequencies different from 1.4 GHz; moreover, the corre¬ 
sponding S '!.4 values are quite low, at the level of a few 
mJy. The optical spectra of all these 10 sources show 
evident and broad emission lines typical of quasars, for 
which a flat radio spectrum is required (see Section [T]). 
In these conditions the impossibility of measuring the 
slope of the radio spectrum prevented them from being 
safely classified as blazars, despite any other require¬ 
ment was fulfilled. 

Blazars can be distinguished from other kind of 
sources through the comparison of the continuum emis¬ 
sion from the active nucleus with respect to the host 
galaxy contribution. A spectroscopic analysis in the op¬ 
tical band is able to show the mutual strength of these 
two compo nents. At this purpose the Ca HfcK br eak 
absorption ( Stocke et al. 199 It Marcha et al. 1996ll . if 
present, is one of the most frequently investigated fea¬ 
tures. A weak non-thermal continuum emission, re¬ 
flecting in a very sharp Ca H&K break, is the foot¬ 
print of a weak central engine and probably shares little 
with the blazar nature, char acterized by much mo re ex¬ 
treme properties. Recently, iMassaro et al.l (|2ni2l l have 
pointed out that, for objects with z < 0.5, the amount 
of the Ca H&K break absorption in the SDSS spectra 
can be determined by computing their {u — r) color 
index from the magnitudes of their SDSS photometry 
adopting the formula: 


{u-r) = {u- r)obs - O.SIA^ 


( 1 ) 


where is the extinction in the r band. In their anal¬ 
ysis, dominant nuclear activity in extragalactic sources 
has been found to correspond to {u — r) values lower 
than a threshold (u — r)* = 1.4. Among the 15 remain¬ 
ing objects only two are at redshift 2 ; <0.5, and their 
(u — r) color index is lower than (u — r)*. 

The list of 15 firmly established blazars selected at 
the end of our investigation is shown in Table [T] Each 
source has been addressed using its SDSS name (col¬ 
umn 1). The angular separation dox between the 
SDSS and ISXPS coordinates is reported (column 2) 
together with the 90% error radius Rx of the ISXPS 



X (A) 


Fig. 1 The SDSS-DR9 optical spectrum of 
SDSS J104031.62-h061721.7. 


source (column 3) provided bv lEvans et al.l (|2014ll . We 
have also reported the angular separation of the SDSS 
source from the PIRST {dpR, column 4) and the WISE 
{doi, column 5) associated counterparts, respectively. 
The values of the radio flux density S' 1.4 (column 6 ) 
have been taken from the most recent version ( 2012 , 
February 16) of the FIRST survey catalog. The val¬ 
ues Sb (column 7) come fro m the Green Bank 6 cm 
Radio Source Catalog ('GBO. fCregorv et"^ 1996 ) with 
the only exception of SDSS J003931.58—111102.4 for 
which the value found in the Parkes-MIT-NRAO sur¬ 
vey fPMN. IWright et al. 1994 ) was reported. The val¬ 
ues of the radio spectral index ar between these two 
frequencies (column 8 ) have been calculated using the 
formula: 


ar = —1.81 • Log (Ss/Sia)- 


( 2 ) 


The values of both the magnitude in the r filter (col¬ 
umn 9) and the redshift 2 : (column 10) come from the 
SDSS-DR9. One of the redshift values, correspond¬ 
ing to SDSS J104031.62-1-061721.7, is reported with a 
question mark as we found a flag in the optical spec¬ 
trum. The flag is due to the fact that the z value 
provided in the SDSS is based on a single, faint emis¬ 
sion line. Assuming that this feature in the optical 
spectrum (Fig. [T]) can be attributed to Mg H we car¬ 
ried out a further check on the redshift measurement. 
The value that we computed, which is reported in Ta¬ 
ble [TJ is very similar to the one (z=0.735) given in the 
SDSS. The values of the unabsorbed flux Fx over the 
0.3-10 keV band (column 11) come from lEvans et all 
(|2014l) . The corresponding values of the X-ray lumi¬ 
nosity Lx (column 12 ) were computed from the for¬ 
mula Lx = i n (dL)"^ ■ Fx, where the luminosity 
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distance was d erived fr o m th e redshift z using the 
web calculatoiH by Wright ( 2006 ). The values that we 
obtained are well above the threshold established 
in the Roma-BZCAT (Section [J), by at least an order 
of magnitude. The quality flag of all the corresponding 
ISXPS sources is equal to 0 meaning a good quality 
X-ray detection (Section [5]). Finally, in column 13 we 
reported the proposed classification as a BL Lac or as 
a FSRQ, mainly driven by the features shown in their 
optical spectrum. 
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Fig. 2 The SED of SDSS J104031.62-h061721.7. Red cir¬ 
cles correspond to radio data at 1.4 GHz and 5 GHz; or- 
ange circles correspon d to infrared da ta from WISE, 2MAS S 
llSkrutskie et al.l[2006l l and UKIDSS llLawrence et al.l[2007ll : 
SDSS data (yellow) are superposed to various data from 
the NED database (brown); X-ray data come from the 
Swift (light green) and the XMM-Newton (green) satellites; 
purple circles correspond to Fermi data. 


4 Discussion 

Blazars are certainly not the only sources characterized 
by a multifrequency emission. In recent years their pe¬ 
culiar properties have been well defined by restrictive 
criteria: the need for establishing them has become 
quite urgent as a consequence of the rapid increase in 
the lists of candidates, a large number of which is sim¬ 
ply the result of automatic cross-matches of the var¬ 
ious all-sky surveys obtained by spacecraft telescopes 
at different energy bands. For this reason, a careful 
check of the overall properties is highly recommended 
to support any candidate before validating the blazar 


^http: //www. astro. ucla.edu/'^wright/CosmoCalc.html 


classification. At the end of our screening procedure, 
that includes some steps that could not be carried out 
automatically, we obtained a list of 15 blazars. We 
verified that the SED of all our blazars is consistent 
with the characteristic double-bumped shape, the for¬ 
mer attributed to synchrotron radiation and the latter 
to inverse Compton emission. One of them is reported 
as example in Fig. [5] All the values of the radio spec¬ 
tral index ar are strictly consistent with a flat radio 
spectrum. According to the classification that we pro¬ 
pose all the sources of our list, with two exceptions, are 
FSRQs with a mean redshift value z = 0.851. FSRQs 
are generally foun d among low synchrotron peak (LSP, 


Abdo et al. 2010l) sources: in these conditions, the de¬ 


tection in the soft X rays corresponds to the inverse 
Compton emission. 



Fig. 3 WISE 3.4 /im sky map centered at the position 
of 2FGL J1040.7-1-0614. The BL Lac object found by our 
procedure, SDSS J104031.62+061721.7, is within the error re¬ 
gion of the 7 -ray source (red crossed line) at variance with 
BZQ J1041+0610, the counterpart associated in the 2LAC. 


In the optical spectrum of SDSS J104031.62+061721.7 
(Fig. [T]) there is evidence of continuum emission with 
no broad emission lines, suggesting a BL Lac classi¬ 
fication. Consid erable variability i n the optical band 
was reported bv iBauer et al.l (|2009l) . We note (Fig. [3]) 
that this source is within the 95% error circle @95 = 
5.88' of the 7 -ray source 2FGL J1040.7+0614; the angular 
separation between the SDSS and the 2FGL coordi¬ 
nates is 3.96'. Conversely, the counterpart 4C -1-06.41 
(named BZQ J1041+0610 in the i?077ia-BZCAT) currently 
associated to this 7 -ray source in the Second LAT 
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AGN Catalog (2LAC, lAckermann et al' 2011) is out¬ 
side of the 2FGL error circle, at a distance of 9.45'; 
for this reason, their association is given with a low 
bayesian probability. Therefore, we consider the even¬ 
tuality that a contribution to the 7 -ray emission cor¬ 
responding to 2FGL J1040.7-1-0614 might be attributed 
to SDSS J104031.62+061721.7, possibly in addition to 
BZQ J1041+0610. 

The other source that does not show broad emission 
lines in its optical spectrum is SDSS J161541.21+471111.7. 
The value that we found for its color index (u — 
r) = 1.34, e ven if lower than the threshold {u — r)* es¬ 
tablished bv iMassaro et al.l (j2012ll . suggests a relatively 
weak non-thermal contribution with respect to the host 
galaxy. Moreo v er, thi s source was recently addressed by 


Massaro et al.l (|2013ll as a possible counterpart of the 


unidentified 7 -ray source (UGS) 2 FGL J1614.8+4703 de¬ 
tected by FerTTii-LAT. 

If we consider the totality of blazars selected by our 
method, including the 152 already reported in the 4*^ 
version of the Aoma-BZCAT, in our list of 167 ob¬ 
jects we found 72 BZBs (~43.1%), 85 BZQs (-50.9%) 
and 10 BZUs (—6.0%), respectively. These rates are 
in agreement with those corresponding to the whole 
blazar population as reported in the Aoma-BZCAT 
(Section [T|). 
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